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ABSTRACT 



Context. Observations of low-mass companions for which the dynamical masses are well constrained help to improve the calibration 
of evolutionary models. Such observations thereby provide more confidence in the estimation of the mass of a companion using the 
photometric methods expected for the next generation of planet finder instruments. 

Aims. The commissioning of a new coronagraph at the Very Large Telescope (VLT) was the occasion to test the performance of 
this technique on the well-known object AB Dor A and its 0.09Mo companion AB Dor C. The purpose of this paper is to refine the 
photometric analysis on this object and to provide an accurate photometric error budget. 

Methods. In addition to coronagraphy, we calibrated the residual stellar halo with a reference star. We used standard techniques for 
photometric extraction. 

Results. The companion AB Dor C is easily detected at 0. 185'' from the primary star, and its magnitudes in H and Ks are in agreement 
with an M5.5 object, as already known from spectroscopic observations. However, these new measurements make the earlier /-band 
photometry less reliable. Finally, the comparison with evolutionary models supports an age of (75 + 25) Myr, contrary to previous 
analyses. These observations demonstrate that coronagraphic observations can be more efficient than direct imaging, not only to 
improve contrast, but also to provide a better photometric estimation as long as a good calibration of the stellar halo is achieved. 

Key words. Stars: individual: AB Dor - Stars: low-mass, brown dwarfs - Techniques: high angular resolution - Methods: observa- 
tional 



1. Introduction 

Since the last decade, the improvement of high angular resolu- 
tion on large telescopes has made possible the discovery of faint 
companions with masses close to the planetary mass regime, 
2M 1207 B (Chauvin et al.ll2Q0 5a): DHTau B dltoh et al. 2005); 
GQLu p B (.NeuhauseTeta^ A B^ic B lchauvin et al. 

\ l2005bh : and CHXR73 B (iLuhman et al.ll2006h . The precise de- 
termination of the mass of most of these companions is presently 
not possible as long as dynamical measurements are missing. 
However, an estimation of the mass can be obtai ned from pho- 
tometric measurements via evolutionary models (iBurrows et al.l 
ll997l:IChabri'ir et al.ll20QQl) . Calibration of these models on very 
low-mass objects, for which the mass is known from other tech- 
niques (radial velocity, astrometry), is highly desi rable to pre- 
pare future instr uments. The inst ruments SPHERE (iBeuzit et al.l 
[2006 ) and GPI ( Macintosh et al.l l2006) will precisely use broad- 
band differential imaging or low-resolution spectroscopy to 
carry out statistical analysis on extrasolar planets. Hence, accu- 
rate mass estimation is critical in this context and a strong effort 
has been made to refine atmospheric models of giant planets and 
brown dwarfs. 

Young, nearby associations are well suited to the identifica- 
tion and the follow-up observations of young dynamical mass 
calibr ators, such as the tight binaries HP 98800 (.Boden et al.l 
I2005h and TWASAab (iKonopackv et al.l l2007h of the TW 
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Hydrae association. Very recently, much attention has been paid 
to the hierachical quadruple system AB Dor, a memb er of 
the eponymous comoving group identified by Zuckerman e t al.l 
(i2004|) . The brightest component AB Dor A was first known as 
a variable star featuring variation of 0. 09 mag in the V b and and 
flares of 0.05 mag near its maximum (llnnis et al. l [T98l . It was 
then recognized as a rapidly-rotating spotted star and was inten- 
sively studied as such in the 1980s. Accurate parallax obtained 
with Hipparcos (tt = (6.92 ± 0. 54) mas, d = (14 .9 ± 0.12) pc; 
iPerrvman et~aD (Il997h ) allowed I Wichmann et af] (11998 ) to de- 
rive a spectral type of Kl, while it was previously thought to be 

a post-T Tauri star. 

A t 9.0'' North the p hysical companion ABDor B (iLimI 

119931: iGuirado et al.ll2QQ6h is r esolved as a tight (A = 0.070 '0 
binarv bv IClose et al.l ('2005'). The object ABDor C is the 
fourth component of this young quadruple system, discovered 
thanks to the reflex motion induced on AB Dor A detected with 
Very Long Baseline I nterferometry and H ipparcos observations 
dOuirado et al.lll997h . IClose et al.1 (l2005h have refined the mass 
estimation of AB Dor C to (0.090 ± 0.005) M© (confirmed later 
by Guirado et al. 2006). A first attempt to image this close and 
low-mass companion with ADONIS at the 3.6m telescope of La 
Silla (ESQ) was unsuccess f ul due to the lack of angular resolu- 
tion (iBoccaletti et al.ll200lh . lciose et al.1 (2005) finally resolved 
AB Dor C at 0.156'' from A using VLT/NACO. They measured 
the near IR absolute magnitudes (see Tab. [21) and a spectral type 
M8 ± 1 for this faint companion. The /- and H- brightnesses 
and the deduced effective temperature were found to be incon- 
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Fig. 1. Coronagraphic image of AB Dor obtained in the Ks filter (Left) and the same image subtracted with a reference star (Right). 
The field of view is 2''. North is up, East is left. Arbitrary false colors are intended to enhance the companion visibility. 



sistent with evolutionary models, considering an age estimate 
of 50^^Q My r based on diff'erent youth indicators. In conclusion, 
IClose et ail (i2005) suggested that theoretical models are actually 
un derestimating the n i ass i n the young age and low-mass regime. 

iLuhman & Potted ('2006) revised the JH K photometry uncer- 
tainties and the spectral type estimation of IClose et al.l (l2005l) 
based on the same data set. Using, in addition, a diff'erent age es- 
timate of 75-150 Myr for the AB Dor association ( Luhman et al. 
l2005l) . they concluded that there was currently no disagree- 
ment between mod els and data. Recen t VLT/SINFONI spectra 
in HK-bands enable Iciose et al ] (I2007h to derive a spectral type 
M5.5± 1 , which confirms the conclusions of Luhman & Potted 
(l2006h . This prolific system illustrates the difficulty of test- 
ing evolutionary model predictions without accurate observables 
(effective temperature and luminosity) and robust age estimate. 

In 2006, a proposal to combine the Simultaneous Diff'erential 
Imaging (SDI) mode of NACO (the Nasmyth Adaptive Optics 
System and Near-Infrared Imager and Spectrog raph) with a 4 
Quadrant Phase Mask coronagraph (iRouan et al. 2000) was ap- 
proved by European Southern Observatory (ESO). During the 
commissioning run we collected data on AB Dor A and C (Sect. 
[21) and here, we present the results of our H and Ks photomet- 
ric analysis (Sect. [3]) together with a detailed estimation of error 
bars (Sect. 131). Results are discussed in Sect.[5l 

2. Observations 

We carried out observations as part of a commissioning run 
on February 16th, 2007 at ESO/Paranal. T he AO-assisted near- 
IR camera NAOS-C ONICA named NACO (iRousset et al.ll2003l: 
iLenzen et al.l l2003l) was equipped w ith a new set of two 4 
Qua drant Phase Masks (iRouan et al.l [2000) to replace the old 
one dBoccaletti et al.ll2004l) . These two masks are operating re- 
spectively in the Ks and H bands, the latter being compati- 
ble with the Simultaneous Diff'erential Imager (SDI) provided 
by the University of Arizona and the Max Planck Institute of 
Heidelberg (Lenzen et al. 2004). 

We observed ABDor (V=6.93, H=4.845, K=4.686, 
Sp=KlIII) with the 4 Quadrant Phase Mask (4QPM) in two 
filters. We obtained 600s (DIT=ls, NDIT=100, Ncycle=6) of 
data in the Ks band (A = 2.18yum, AA = 035jim) on the target 
and a similar integration time on a reference star (HD41371) 
chosen with the same visible and IR fluxes (V=7.10, K=4.724, 
Sp=KOIII) and observed at the same parallactic angle. This 



optimal observing strategy preserves the orientation of the 
telescope pupil with respect to NACO, and therefore, reduces 
the difl'erential aberrations between the star and its reference as 
the matching of spider spikes in the two images. In addition, we 
observed AB Dor with SDI (the 4QPM being installed in the 
beam) for 936s (DIT=8s, NDIT=13, Ncycle=9). Instead of a 
reference star, the calibration of the speckled halo is obtained 
simultaneously in difl'erent filters (A = 1.575, 1.600, 1.625yum, 
AA = 0.025yum). However, a second level of calibration is 
required to reduce the impact of diff'erential aberrations and 
the speckle chroma ticity (as the phase varies with wavelength, 
iMarois et al. For this purpose, we obtained two ob- 

servations with the field of view rotated by 60°. Diff'erential 
aberrations are assumed static in this case. 

Coronagraphic observations with NACO are preceded with 
an acquisition template that provides an out-of-mask PSF to be 
used as a photometric reference. However, this reference is ob- 
tained with a diff'erent setup than coronagraphic frames. Because 
AB Dor is a bright object, a Neutral Density (ND) is needed to 
avoid detector saturation (for Ks band data only), and as default 
a full aperture stop (Full) is used instead of a stopped aperture 
(FulLuszd) as for coronagraphic templates. Photometric mea- 
surements have to be corrected from these values. 

Seeing conditions and AO correction were good. The aver- 
age seeing was (0.91 ± 0.13)" for the Ks observing block and 
(0.78 ± 0.12)" for the SDI data while the coherent energy mea- 
sured from residual slopes of the AO system was, respectively, 
(46 ± 5)% and (52 ± 4)%. 



3. Data reduction and photometric measurements 

We processed the data with standard reduction routines to correct 
for bad pixels, flat field uniformity and to subtract an average sky 
background. We obtained flat fields without the coronagraph in 
the beam, although it would have been better to reduce the pres- 
ence of dust particles on the substrate. However, the field of view 
slightly drifted in front of the detector plane (as the instrument 
rotates at the Nasmyth focus), making the precise registration of 
dust features on the coronagraph substrate impossible. 

We corrected for the ND attenuation on the Ks data (a fac- 
tor of 89) and for the diff'erence of pupil stop (a factor of 0.808 
estimated from geometrical assumptions) to build a master nor- 
malized PSF. The individual coronagraphic cycles were then co- 
added for the star and the reference separately so as to provide 
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Fig. 2. Coronagraphic image of AB Dor obtained with SDI (Left) and the same image subtracted with a field-rotated AB Dor image 
(Middle) compared to the 2 wavelengths subtraction (Right). Subtracted frames show a positive image of the companion (yellow) 
and a negative one (red). The SDI frame is not exploitable because of the small angular separation making the positive and negative 
component self- subtracting. The field of view is 2'\ North is up, East is left. Arbitrary false colors are intended to enhance the 
companion visibility. 



two images. Recentering was not required. At this stage, the 
companion is yet visible (Fig. [T] and Fig.O Left), but we needed 
the subtraction of the reference star to remove the PSF halo and 
to accurately measure the photometry. To estimate the intensity 
factor between the coronagraphic image of the star and that of 
the reference, we considered several methods: total intensity ra- 
tio, balance of positive and negative fluxes, balance of positive 
and negative pixels, and minimization of the total residual in- 
tensity. We adopt the last method as a baseline paying attention 
to the presence of the companion to avoid a bias. We found an 
intensity factor of 1.04 ± 0.01 with small dispersion between 
methods and produced a subtracted coronagraphic image (Fig. 
[D Right). 

Once the stellar contribution is removed at the companion 
location a thorough estimation of the companion intensity be- 
comes possible. Here also, we compared several methods like 
aperture photometry, PSF fitting with a 2D gaussian, PSF to 
companion maximum intensity ratio, and minimization of the 
residual after PSF subtraction (on the companion). In the two 
first cases, the companion intensity is integrated in a limited 
aperture (a few pixels in radius) and a correction of -O.lmag is 
needed to account for the intensity in the PSF wings (comparison 
of encircled energy). As a baseline, we used the minimization to 
derive the average photometry. 

In the particular case of SDI data, two wavelengths subtrac- 
tion is inappropriate since AB Dor C does not contain methane 
and the bifurcation point (the distance at which the companion 
image in the res caled frame falls at diff'erent pixels than in the un- 
rescaled frame. Pftiatte et al.ll2007h is located at 3.2''. However, 
observational procedure with SDI requires the acquisition of a 
field-rotated image of the same star to correct for the diff'erential 
static aberrations and the chromatic aberrations inherent to the 
SDI technique. Therefore, we measured H band flux of AB Dor 
C with the field-rotated image as a reference star both in a single 
SDI filter (A A = 0. 025 fim) and with the combination of all filters 
(AA = 0.075yum) to approach the broadband magnitude. Hence, 
two values are provided for the H band magnitudes that we note: 

As a result, we found the following magnitude diff'erences : 
AntHu =4.71, ArriH^^ = 4.62 and ArriKs = 4.56. The angular sep- 
aration is 0. 1 85 consistent with the orbital solution presented in 



Table 1. Photometric uncertainties in magnitude. 



sources of error 


H 


Ks 


pupil stop 


0.05 


0.05 


ND 




0.06 


intensity factor 


0.05 


0.04 


method for extraction 


0.13 


0.13 


aperture size 


0.04 


0.01 


distance 


0.03 


0.03 


2MASS 


0.03 


0.02 


filters conversion 


0.03 


0.01 


total quadratic error 


0.16 


0.16 



iNielsen et al.l (|2005|) . However, it was not the goal of this paper 
to discuss the astrometry of the companion since it was already 
characterized from Hipparcos data and confirmed by the afore- 
said papers. In addition, the presence of the coronagraph makes 
the estimation of the astro metry less accurate unless appropriate 
techniques are considered (iMarois et al.ll2006l) . 

4. Photometric errors 

We identified several sources of errors in the photometric extrac- 
tion that we analyze in this section to derive error bars: 

- The pupil Lyot stop correction made on the PSF flux. As 
mentioned, a geometrical comparison of the Full aperture 
and the FulLuszd aperture (undersized by 10%) leads to a 
correction of 0.808. However, a photometric measurement 
obtained in a previous observing run suggests a value of 
0.775. Then, an uncertainty of 4% was considered. 

- The ND correction made on the PSF flux (only for data in 
Ks). Our measurement in the Ks band yields an attenuation 
factor of 89 + 3.6. 

- Intensity factor between the coronagraphic image of the star 
and that of the reference. The variety of method we used to 
measure this parameter is providing a good estimate of the 
uncertainty. The precision achieved is 1%. 

- Photometric extraction. Here again, we used a variety of 
methods to consolidate the result. However, we identified it 
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Table 2. Photometry of AB Dor, A being the primary and C the companion converted into the 2MASS system. 



UT Date 


Mj 
(mag) 


(mag) 


(mag) 


A 

(mas) 


SpT 


(K) 


References 


04.02.2004 

07.01.2005 
24.01.2006 
16.02.2007 
16.02.2007 


^•"^-0.24 


9 17+0.13 

Q Ql+Sil 
^"^^-0.31 

8.71+0.16 
8.64 + 0.16 


o <O+0.12 

8 92+^-28 
8.63 + 0.17 

8.38 + 0.16 


156 
156 
220 
200 
185 


M8 + 1 
M6 + 1 

M5.5 + 1 


2600:|| 
2840^1^^ 

2925+170 


Close et al. (2005) 
Luhman & Potter (2006) 
Close et al. (2007) 

Thatte et al. (2007). Close et al. (2007) 
this work 
this work 



as the major source of uncertainty in our measurement. The 
dispersion between methods is of about O.lSmag. 

- Aperture size, when aperture photometry is used for pho- 
tometric extraction. The radius of the aperture photometric 
mask is set to 1. 22 A /Diyot, Diyot being the diameter of the 
Lyot stop instead of that of the telescope (7.2m instead of 
8m). A variation of 1 pixel on this radius provides the error 
bar. Although aperture photometry is not used to extract the 
photometry (see Sect. [3]) it actually enters in the error term 
"method of extraction", but has a wavelength dependence ac- 
cording to the PSF sampling. 

- The star distance which is known to an accuracy of 0.12 pc. 

- The 2MASS photometric uncertainty for the primary star 
magnitude. 

- The conversion between systems of magnitude (NACO to 
2MASS to GIT). 

Table [T] gives the correspondence of these uncertainties in 
magnitude. To derive the AB Dor C absolute magnitudes in the 
H and Ks 2MASS system, we have estimated the NACO to 
2MASS filters transformations based on the primary and sec- 
ondary spectral types and the filter transmission curves. The 
transformations found are (-0.02+0. 03)mag, (-0.04±0.03)mag, 
and (0.00±0.01)mag in H^, Hi,a, and Ks respectively. Assuming 
errors are added quadratically and a distance of J = (14.9 ± 
0.12) pc, the new H- and Ks- absolute magnitudes of AB Dor 
C were obtained with uncertainties and are given with the results 
of previous works in Tab. |2l 




0.10 1.00 
Age (Gyr) 



Fig. 3. Interpolated evolutionary tracks of "Chabrier et alj (|200Q|) 
for 0.085 (dashed), 0.090 (solid) and 0.095 (dotted) solar masses 
in the range 0.01-10 Gyr comp ared with our phot ometric mea- 
surements (red box) those of iGlose et aP (l2QQ5l) (green dia- 
mond) and IClose et al.l (l2QQ7l) (purple triangle) and those of 
iLuhman & Potted (|2006) (blue circle). Upper and lower plots 
are respectively for Ks and H photometry (in the CIT sys- 
tem of magnitudes). The left plots are for the first age scenario 
(75 ± 25 Myr) and the right plots for the second (120 + 20 Myr). 



5. Discussion 

Our new H and Ks measurements of ABDor C and the de- 
rived H - Ks = 0.29 + 0.23 color are consistent with the recent 
M5.5 + 1 spectral type estimation of IClose et al.l (l2007l) (using 
spectroscopic analysis rather than photometry). These new mea- 
su rements also weaken th e reliability of the /-band photometry 
of IClose et aP (l2QQ5l) and iLuhman & Potted (l2QQ6b obtained in 
poor a tmospheric condit ions and that "may be systematically too 
faint" (IClose et al ]l2007h . The J - Ks would be indeed even red- 
der than before and surprisingly red compared to M5.5 objects 
of the Pleiades. Although not impossible, the existence of an IR 
excess for AB Dor C seems to be unlikely as AB Dor A (Kl), 
Ba (M3.5 + 1.5), and Bb (M4.5 + 1.5) do not present abnormal 
(/ - Ks) colors. 

Our new H and Ks absolute magnitudes of AB Dor C (color 
corrected again between the 2MASS and CIT systems) are 
compared with previous measure ments and ov e r-plott ed on 
evolutionary model predictions of Chabrier et al.' (2000) for a 
dynamical mass of (0.09 + 0.005)Mo (Guirado et al. 2006). For 
the following analysis, we have considered two scenarii of ages 
for the AB Dor system, (75 + 25) Myr and (120 + 20) Myr, as 



the determination of this fundamental parameter is still debated. 

The first age estim ate of ~ 50 Myr was given by 
IZuckerman et al.l (l2004l) for the whole AB Dor association. We 
used multiple youth indicators, such as Ha emission, strong 
lithium 6708 A absorption, large v sin/, large X-ray flux, and a 
location above the main sequence of a color (V - K) magni- 
tude (M^) diagram for th r ee mid -M type members. IClose et aP 
(2005) and ' Nielsen elaP (l2Q05l) estimated a more precise age 
range to 50^2o based on comp arisons with other yo ung clus- 
ters properties. Similar analysis led lLopez-Santiago et a l. (2006) 
to confirm this age estimate. More recently. Uanson et al. (2007^) 
focused their age-dating criteria on the the two M dwarfs AB Dor 
Ba and Bb, comparing their eff'ective temperatures to evolution- 
ary track predictions, which are relatively robust for low-mass 
stars. However, uncorrect spectral types for both components led 
them to over constrain the upper age limit for this system. We 
show model predictions on Fig. [3] (Left), and compared to our 
measurement s for an age o f (75 + 25) Myr. Contrary to the mea- 
surements of Close et^.l ((2005) and Luhman & Potter (2006) 
obtained in less favorable conditions (smaller angular separation 
and lower AO correction), our results are clearly in good agree- 
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ment with evolutionary tracks in both H and Ks and would con- 
firm the validity of evolutionary models in this range of masses 
and ages. 

Us ing comparable age-dating indicators, iLuhman et alJ 
(l2QQ5h derived a second estimate of 75-150 Mvr for the AB Dor 
association. Moreover, like llnnis et aO (Il986h . they have sug- 
gested that the AB Dor association might belong to the Pleiades 
supercluster sharing common kinematic origin. Based on the 
3D dynamical e volution analysis o f the ABDor association 
and the Pleiades, lOrtega et"an (1200 7^ obtained a similar age of 
(120 + 20) Myr consistent with both groups being coeval. If we 
now consider this age- scenario shown on Fig. [3] (Right), our 
measurements are not well fitted in both bands by the model 
predictions. The model predictions actually underestimate the 
observed luminosity of AB Dor C by a factor of ~ 1.6 (0.5 mag), 
and thus overesti mate its mass (opposite of what was found by 
IClose et al ] l2005h . 

Future age characterization and confirmation of our photo- 
metric results should help to draw a robust conclusion to esti- 
mate the accuracy of the model predictions in this range of ages 
and masses. More interestingly, the results presented here evi- 
dence the ability of broadband photometry combined with coro- 
nagraphy to retrieve the actual companion mass, providing opti- 
cal quality is met to allow the use of a coronagraph. This conclu- 
sion is important in the context of planet finder instruments like 
SPHERE and GPI, using diff'erential imaging observations of 
some planetary spectral features and for which planetary masses 
will be assessed via the comparison of spectral contrasts and evo- 
lutionary models. 

Acknowledgements. We would like to thanks the ESO team at the telescope who 
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